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£?■ Abstract 



An interpretation of the formation of halo in accelerators based on quantum- 



CN \ like theory by a diffraction model is given in terms of the transversal beam 



motion. Physical implications of the longitudinal dynamics are also examined. 
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'^' I. INTRODUCTION 

Q-c Recently the description of the dynamical evolution of high density beams by using the 

collective models, has become more and more popular. A way of developing this point of 
K^ ! view is the quantum-like approach [0] where one considers a time-dependent Schrodinger 

H \ equation, in both the usual linear and the less usual nonlinear form, as a fluid equation 

for the whole beam. In this case the squared modulus of the wave function (named beam 
wave function) gives the distribution function of the particles in space at a certain time 0. 
The Schrodinger equation may be taken in one or more spatial dimensions according to the 
particular physical problem; furthermore the motion of the particles in the configuration 
space can be considered as a Madelung fiuid if one chooses the equation in its linear version. 
Although the validity of the model relies only on experiments and on the new predictions 
which must be verified experimentally, we like to invoke here a theoretical argument that 
could justify the Schrodinger quantum-like approach. Let us think of particles in motion 
within a bunch in such a way that the single particle moves under an average force field 

1 



due to the presence of all others and collides with the neighbouring ones in a complicated 
manner. It is obviously impossible to follow and describe all the forces deterministically. One 
then faces a situation where the classical motion determined by the force-field is perturbed 
continuously by a random term, and one finds immediately a connection with a stochastic 
process. If one assumes that the process is Markovian and Brownian, one easily arrives at 
a modification of the equations of motion in such a manner that would be synthesized by 
a linear Schrodinger equation depending on a physical parameter that has the dimension of 
action PQ . Wave quantum mechanics follows if this parameter coincides with the Planck's 
constant h, whereas the quantum-like theory of beams is obtained if one chooses it as the 
normalized emittance e |1]] . In both cases, the evolution of the system is expressed in terms 
of a continuous field ip which defines the so-called Madelung fluid. We may notice that the 
normalized emittance e with the dimension of an action is the natural choice in the quantum- 
like theory, that finds the analogue in the Planck's constant h because it reproduces the 
corresponding area in the phase-space of the particle. 

We here point out that, after linearizing the Schrodinger-like equation, for beams in 
an accelerator, one can use the whole apparatus of quantum mechanics, keeping in mind 
a new interpretation of the basic parameters (for instance the Planck's constant h — > e 
where e is the normalized beam emittance). In particular one introduces the propagator 
K {xf,tf\xi,ti) of the Feynman theory for both longitudinal and transversal motion. A 
procedure of this sort seems effective for a global description of several phenomena such as 
intrabeam scattering, space-charge, particle focusing, that cannot be treated easily in detail 
by "classical mechanics". One consequence of this procedure is to obtain information on 
the creation of the Halo around the main beam line by the losses of particles due to the 
transversal collective motion. 



II. TRANSVERSAL MOTION 

Let us indeed consider the Schrodinger like equation for the beam wave function 

iedti; = -^dl^P + U{x,t)^P (1) 

in the linearized case U {x,t) does not depend on the density |^| . e here is the normalized 
transversal beam emittance defined as follows: 

e = moC7/5e , (2) 

e being the emittance usually considered, (we may also introduce the analogue of the de 
Broglie wavelength as A = e/p). Let us now focus our attention on the one dimensional 
transversal motion along the a;-axis of the beam particles belonging to a single bunch and 
assume a Gaussian transversal profile for particles injected into a circular machine. We want 
to try a description of interactions that cannot be treated in detail, as a diffraction through 
a slit that becomes a phenomenological boundary in each segment of the particle trajectory. 
This condition should be applied to both beam wave function and beam propagator K. The 
result is a multiple integral that determines the actual propagator between the initial and 
final states in terms of the space-time intervals due to the intermediate segments. 



r+b 

K {x + xo,T + t\x', 0) = / K {x + Xq, t\xo + yn,T + {n - l)r') 

J ~b 

xK{x + yn,T+{n- l)r'|xo + Vn-uT + (n - 2)r') 

X ■■ ■K{x + yi,T\x', 0) dyidy2 ■ ■ ■ dyn (3) 

where r = nr' is the total time spent by the beam in the accelerator (total time of revolutions 
in circular machines), T is the time necessary to insert the bunch (practically the time 
between two successive bunches) and (—6, +h) the space interval defining the boundary 
mentioned above. Obviously h and T are phenomenological parameters which vary from a 
machine to another and must also have a strict correction with the geometry of the vacuum 
tube where the particles circulate. 

We may consider the two simplest possible approximations for K {n\n — 1) = 
K (xo + yn,T +{n- 1)t'\xq + |/„„i + {n - 2)r'): 

1. We substitute the correct K with the free particle Kq assuming that in the r' interval 
(r' <^ r) the motion is practically a free particle motion between the boundaries 

(-6, +6). 

2. We substitute it with the harmonic oscillator K^^ {n\n — 1) considering the betatron 
and the synchrotron oscillations with frequency ti;/27r 

III. FREE PARTICLE CASE 

We may notice that the convolution property (|^) of the Feynman propagator allows us 
to substitute the multiple integral (that becomes a functional integral for n — > cxd and 
t' — > 0) with the single integral 

r+b 

i^ (x + Xo,T + r|x', 0) = / dyK {x + Xq,T + t\xo + y,T)K {xo + y,T\x' ,Q)dy (4) 

J-h 



After introducing the Gaussian sht exp 

K {x + Xq,T + r|x', 0) 
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instead of the segment (—6, +h) we have 
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where Vq = ^°~^ and XqIs the initial central point of the beam at injection and can be 
chosen as the origin {xq = 0) of the transverse motion of the reference trajectory in the 
frame of the particle, h must be interpreted as the normalized beam emittance in 
the quantum-like approach. 

With an initial Gaussian profile (at t = 0), the beam wave function (normalized to 1) is 
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- being the r.m.s transversal spot size of the beam; the final beam wave function is: 
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The final local distribution of the beam that undergoes the diffraction is therefore 

p{x) = |0(a;)| = BB*exp —ax"^ 
where a = —{C + C*) and the total probability per particle is given by 

P= I dxp{x) =BB\ 



(9) 



(10) 



Under certain physical conditions (such as the LHC transversal, Table-I), P ~ ^^. 



IV. OSCILLATOR CASE 

Similarly we may consider the harmonic oscillator case (betatronic oscillations and syn- 
chrotronic oscillations) to compute the diffraction probability of the single particle from the 
beam wave function and evaluate the probability of beam losses per particle. The propagator 
K^ (x, T + r||/, T) in the later case is: 
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where 
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Under some physical situations (such as the LHC transversal case) we have, Puj ~ 
'Ta'^ sintuiT) - ^"^ *^^ approximate formulae for P and P^), when applicable, the parameter r 
does not play a significant role. 



V. LONGITUDINAL MOTION 

As far as the longitudinal motion is concerned the quantum-like approach appears to 
be quite appropriate to obtain information on the modified length (and consequently the 
stability) of the bunches both in the linear and circular accelerators.. To be more specific it 
describes a large number of important nonlinear phenomena that are present in RF particle 
accelerators (with residual addition of longitudinal coupling impedance) as well as in cold 
plasmas 0. 

We introduce the Gaussian parameter 6, as we did with the Gaussian slit e~^ /'^^ in the 
transversal motion and look for a phenomenological solution of the equation for the beam 
wave function ifj 
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-dl^ + -mouj^x^^ + A 
2y^mQ 2 



(17) 



where cu is the synchrotron frequency, A represents the coupling with non-linear terms and 
X is the longitudinal particle displacement with respect to the synchrotronous one. 

The Feynman propagator is given by Eq. (|ll|) and the initial wave function can be again 
assumed as a Gaussian wave packet. The main difference with the transversal case stays in 
the numerical values of the parameters that exhibit a different physical situation and require 
a different physical interpretation. 



VI. PRELIMINARY ESTIMATES 



Examples of the numerical calculations for two projects (LHC for ions and HIDIF for 
heavy ions) with very different physical characteristics are reproduced in the following tables. 

TABLE-I: Circular Machines: Transversal Case 

i) HIDIF (storage ring) 



Parameters 




LHC (at injec 


Normalized Transverse Emittance 


3.75 mm mrad 


Total Energy, 


E 


450 GcV 


1 




1.2 mm 


T 




25 nano sec. 


T 




88 sec. 


h 




1.2 mm 


1 




1.41 X 10^ m 


p 




3.39 X 10-5 


UJ 




4.44 X 10^ Hz 


1 




1.03 X 10^ m 


P. 




3.40 X 10-5 



13.5 


mm ] 


iirad 


5Gev 




1.0 mm 




100] 


nano 


sec. 


4.66 


sec. 




1.0 mm 




1.96 


X 10^ 


m 


2.37 


X 10- 


-3 


1.15 


X 10^ 


Hz 


2.07 


X 10- 


-im 


3.00 


X 10- 


-3 



TABLE-H: Circular Machines: Longitudinal Case 



Parameters 

Normalized Longitudinal Emittance 

Total Energy, E 
1 

\/OL 

T 

T 
b 

UJ 



P... 



LHC (at injection) 

1.00 eV sec. 

450 GeV 

7.7 cm 

25 nano sec. 

88 sec. 

7.7 m 

4.23 X 10^ Hz 

1.14 X 10^ m 

0.575 



TABLE-HI: RF Main LINAC of HIDIF 



Parameters 

Normalized Longitudinal Emittance 
Total Final Energy, E 



T 

T 
b 

UJ 



P., 



0.7 keV nano sec. 

5 Gev 

15 cm 

75 micro sec. 

4.9 X 10-^ sec. 

15 m 

4.13 X 10^ Hz 

6.72 X 10-2 m 

0.707 



The machine parameters of tables I, II and III are derived from 0, 0. In particular uj of 
Table-Ill is calculated on the basis of the "Main LIN AC" Table (page 198 of 0) with the 
standard formula: 



2 eEujRp sin 



^ = ^Z^^, (18) 

where the symbols have the usual meaning. 



VII. COMMENTS AND CONCLUSIONS 

Transversal Motion: This use of a quantum-like approach appears a simple powerful tool 
for the analysis of the evolution of a beam in linear and circular accelerators and storage 
rings. 

Indeed the introduction of a very limited number of phenomenological parameters (in 
our simplified model the only parameter b) in the beam quantum-like equations and the use 
of the Schrodinger-type solutions allow us to calculate how the bunches evolve and modify 
owing to the forces (linear and non-linear) acting on the particles. 

As far as the betatronic oscillations are concerned the mechanism of the diffraction 
through a slit appears a very adequate phenomenological approach. Indeed we can interpret 
the probability (local and total) for a particle leaving its position as the mechanism of 
creating a halo around the main flux. 

The values of r, u are strictly connected with the characteristic parameters of the designs 
of the accelerators (in our example LHC and HIDIF) 

The phenomenological parameter b represents several fundamental processes that are 
present in the beam bunches (and play a determinant role in the creation of the halo) such 
as intrabeam scattering, beamstrahlung, space-charge and imperfections in the magnets of 
the lattice that could cause non-linear perturbative effects. 

We like to recall here the analogy with the diffraction through a slit in optics where it 
represents a much more complicated physical phenomenon based on the scattering of light 
against atomic electrons. 

r is the total time spent in the accelerator by a single bunch, T may coincide with 
the average time interval between two successive injections and u is the betatronic average 
frequency given by 2TTQfr, fr being the revolution frequency. 

The fact that a small number of parameters can take into account many physical processes 
is a very nice feature of the quantum-like diffraction approach. However the deep connection 
between this method and the actual physical process as well as the nonlinear dynamical 
classical theory is necessary to be understood. 

We remark now the following points 

1. The total probability (per particle) calculated from the free particle propagator (P) 
and from the harmonic oscillator one (P^^) appear very near for the two different 
circular systems, LHC and HIDIF. 

2. The local distribution between the two however looks quite different for the free and 
harmonic oscillator case, thus giving us a profile of the halo which appears particularly 
interesting in the HIDIF case (final Gaussian width ~ -7^ ~ 2.07 x 10~^ m) 



3. The HIDIF scenario, as we expect because of the higher intensity, exhibits a total 
loss of particles (and beam power) which is at least 10^ times higher than LHC The 
picture we have obtained for the transversal motion in the two analyzed examples (on 
the basis of the parameters provided by the latest designs) is encouraging because the 
halo losses are under control. In both cases the estimated losses of the beam power 
appear much smaller than the permissible 1 Watt/m. 

Longitudinal motion The formulae (|^ and ( p!B| ) can be used for calculating the motion 
of the length of the bunch related to the synchrotron oscillations in both linear and circular 
machines. In this case we must consider only the propagator of the harmonic oscillator 
which is the simplest linear version of the classical dynamical motion for the two canonical 
conjugate variables that express the deviations of an arbitrary particle from the synchronous 
one namely the RF phase difference A0 = — 0^ and the energy difference AE = E — E^. 
Our examples are again the LHC synchrotron oscillations and the ones of the main LINAC 
in the HIDIF project. The phenomenological Gaussian function e~^ ^^^ acquires a different 
meaning from the one it had in the transversal motion. Our analysis deals with a Gaussian 
longitudinal profile and predicts a coasting beam in LHC and a quite stable bunch in the 
main LINAC of HIDIF. 

We may therefore conclude that our approach although preliminary is interesting and 
particular attention is required in treating the longitudinal motion where the nonlinear 
space-charge forces are very important. So the quantum-like method appears promising for 
the future simulations in beam physics. 
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